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Abstract 

Differential cross sections for dijet photoproduction and this process in associ- 
ation with a leading neutron, e + + p — > e + + jet + jet + X (+n), have been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
40 pb^ 1 . The fraction of dijet events with a leading neutron was studied as a 
function of different jet and event variables. Single- and double-differential cross 
sections are presented as a function of the longitudinal fraction of the proton 
momentum carried by the leading neutron, xl, and of its transverse momentum 
squared, p\. The dijet data are compared to inclusive DIS and photoproduction 
results; they are all consistent with a simple pion-exchange model. The neutron 
yield as a function of xl was found to depend only on the fraction of the proton 
beam energy going into the forward region, independent of the hard process. No 
firm conclusion can be drawn on the presence of rescattering effects. 
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1 Introduction 



The transition of an initial-state proton into a final-state neutron, p — > n, has been 
extensively studied in hadronic reactions [1—7]. A successful phenomenological description 
of these results uses Regge theory and interprets the interactions as an exchange of virtual 
isovector mesons, such as it, p, and 02 [8-11]. At small values of the squared momentum 
transfer, t, between the proton and the neutron, the p — > n transition is expected to be 
dominated by the exchange of the lightest meson, the pion. 

Leading baryon processes have been previously studied in ep collisions at HERA [12-20]. 
Some of these studies were performed involving a hard scale, such as the virtuality of the 
photon exchanged at the lepton vertex, Q 2 , in deep inelastic scattering (DIS) [12,13,16,20]; 
the jet transverse energy, , in photoproduction of dijets [14]; or the charm mass in 
heavy-flavor production [17]. 

Even though a hard scale is involved, the p — > n transitions are still expected to be 
dominated by pion exchange. The cross section for this type of process in ep collisions 
can be written as 

d df = Jtt/pOl, t)a e ^ex(s ). (1) 

This formula expresses the Regge factorization of the cross section into the pion flux 
factor f n / p (xL,t), which describes the splitting of a proton into an n-n system, and the 
cross section for electroproduction on the pion, a en ^ e x(s'). Here, xl is the fraction of the 
incoming proton beam energy carried by the neutron, and s and s' = (1 — xl)s are the 
squared center-of-mass energies of the ep and of the eir systems, respectively. 

Comparisons between neutron-tagged and untagged cross sections provide tests of the 
concept of vertex factorization [21]. Under this hypothesis, the shape of the distribution 
of some photon variable V would neither depend on the presence of a neutron nor explicitly 
on its kinematic variables xl and t. Similarly, the xl and t spectra of the neutrons would 
be independent of the photon variable V. The cross section can then be written as 

^=ff§^ = 9(*L,t)G(V), (2) 

where g(xi,t) and G(V) are functions of the neutron and photon variables respectively. 
The Regge factorization expressed in Eq. ([1]) violates this vertex factorization because a en 
has different s' dependences for different processes and s' depends on xl- This will be 
further explained in Section [71 and violations of vertex factorization are therefore to be 
expected. 

Rescattering effects, where the baryon interacts with the exchanged photon [22-25], are 
expected to increase with increasing size of the virtual photon, i.e. decreasing Q 2 . This 
was observed in a measurement of leading neutrons in DIS and photoproduction [20]. 
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In high-.Ey* jet photoproduction with Q 2 ps 0, two types of processes contribute to the 
cross section, namely direct and resolved photon processes. In direct processes, the ex- 
changed photon participates in the hard scattering as a point-like particle. In resolved 
processes, the photon acts as a source of partons, one of which interacts with a parton 
from the incoming hadron, see Fig. HJ The more complex structure of the resolved pho- 
ton may increase the probability for the leading baryon to rescatter. This can cause the 
baryon to be scattered out of the detector acceptance, resulting in a depletion of detected 
baryons. Thus, fewer leading baryons (i.e. more rescatterings) are expected in resolved 
than in direct processes. 

This effect was searched for, but not confirmed, in diffractive production of dijets in 
photoproduction [26,27] and DIS [26,28], where the leading proton has Xl ~ 1- However, a 
comparison of leading neutron rates in photoproduction and DIS showed a scale dependent 
suppression of neutrons [14,17,20]; the rates of neutrons were in good agreement with the 
expectations from rescattering models [22,23]. 

This paper reports the observation of the photoproduction of dijets in association with a 
leading neutron: 

e + +p -y e + + jet + jet + X + n, (3) 

where X denotes the remainder of the final state. The number of events is almost an order 
of magnitude higher than used for previous results [14, 18]. Cross sections are presented 
as functions of the jet transverse energy, E^ , jet pseudorapidity, if et , the fraction of 
the photon energy carried by the dijet system, i° BS , the photon-proton center-of-mass 
energy, W, and the fraction of the proton four-momentum participating in the reaction, 
x° BS . In addition, the fraction of photoproduction events with a leading neutron as 
functions of these variables is shown as a test of vertex factorization. Finally, the xl and 
p\, distributions of the leading neutrons are shown in dijet photoproduction and compared 
to similar results in DIS [20]. 

2 Experimental setup 

The data sample used in this analysis was collected with the ZEUS detector at HERA and 
corresponds to an integrated luminosity of 40pb _1 taken during the year 2000. During 
this period, HERA operated with protons of energy E p = 920 GeV and positrons of energy 
E e = 27.5 GeV, yielding a center-of-mass energy of yfs = 318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [29]. A brief outline 
of the components most relevant for this analysis is given below. Charged particles were 
tracked in the central tracking detector (CTD) [30], which operated in a magnetic field of 
1.43 T provided by a thin superconducting coil. The CTD consisted of 72 cylindrical drift 
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chamber layers, organized in 9 superlayers covering the polar-angku region 15° < 9 < 164°. 
The transverse-momentum resolution for full-length tracks was ct(pt)/pt — 0.0058pr © 
0.0065 © 0.0014/p r , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [31] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy 
resolutions, as measured under test-beam conditions, are a(E)/E = 0.18/\^E for elec- 
trons and a(E)/E = 0.35 /VE for hadrons (E in GeV). The forward-plug calorime- 
ter (FPC) [32] around the beam-pipe in the FCAL extended calorimetry to the region 
rj m 4.0 — 5.0. It was a lead-scintillator calorimeter with a hadronic energy resolution of 
a(E)/E = 0.65/v 7 ^ © 0.06 (E in GeV). 

The forward neutron detectors are described in detail elsewhere [20,33]; the main points 
are summarized briefly here. The forward neutron calorimeter (FNC) was installed in 
the HERA tunnel at 9 = 0° and at Z = 106 m from the interaction point in the proton- 
beam direction. It was a lead-scintillator calorimeter, segmented vertically into towers 
to allow the separation of electromagnetic and hadronic showers by their energy sharing 
among towers. The energy resolution for neutrons, as measured in a beam test, was 
<j(E n ) / E n = 0.70/v^^, with neutron energy E n in GeV. The energy scale of the FNC was 
determined with a systematic uncertainty of ±2%. The forward neutron tracker (FNT) 
was installed in the FNC at a depth of one interaction length. It was a hodoscope designed 
to measure the position of neutron showers, with two planes of scintillator fingers used 
to reconstruct the X and Y positions of showers. The position resolution was ±0.23 cm. 
Veto counters were used to reject events in which particles had interacted with the inactive 
material in front of the FNC. Magnet apertures limited the FNC acceptance to neutrons 
with production angles less than 0.75 mrad, which corresponds to transverse momenta 
Pt < E n 9 max = 0.69 xl GeV. 

The luminosity was determined from the rate of the bremsstrahlung process, ep — > e'yp, 
where the photon was measured with a lead-scintillator calorimeter [34, 35] located at 
Z = -107 m. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
center of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as r\ — — In (tan f), where the polar angle, 9, is measured with respect to the proton beam 
direction. 
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3 Data selection and kinematic variables 



A three-level trigger system was used to select events online [29,36]. At the second level, 
cuts were made to reject beam-gas interactions and cosmic rays. At the third level, jets 
were reconstructed using the energies and positions of the CAL cells. Events with at least 
two jets with transverse energy in excess of 4.5 GeV and \rf ct \ below 2.5 were accepted. 
No requirement on the FNC was made at any trigger level. 

Offline, tracking and calorimeter information were combined to form energy-flow objects 
(EFOs) [37,38]. The jp center-of-mass energy, W, was reconstructed using the Jacquet- 
Blondel method [39] as Wjb = a/^jbs, where i/jb = ^(-^i — Ez,i)/2E e is an estimator of 
the inelasticity variable y, and Ez,i = Ei cos Of, Ei is the energy of EFO i with polar angle 
9i. The sum runs over all EFOs. The energy Wjb was corrected for energy losses using 
the Monte Carlo (MC) samples described in Section 4. After corrections, the sample was 
restricted to 130 < W < 280 GeV. Events with a reconstructed positron candidate in the 
main detector were rejected. The selected photoproduction sample consisted of events 
from ep interactions with Q 2 < 1 GeV 2 and a mean Q 2 « 10 -3 GeV 2 . 

The kj< cluster algorithm [40] was used in the longitudinally invariant inclusive mode [41] 
to reconstruct jets in the measured hadronic final state from the energy deposits in the 
CAL cells (calorimetric jets). The axis of the jet was defined according to the Snowmass 
convention [42]. The jet search was performed in the (77 — <p) plane of the laboratory 
frame. Corrections [43] to the jet transverse energy, E^ , were applied as a function of 
the jet pseudorapidity, if et , and E^ , and averaged over the jet azimuthal angle. Events 
with at least two jets of E^ 1 ^ > 7.5(6.5) GeV, where E^ is the transverse energy of 
the highest (second highest) E 3 ^ jet, and — 1.5 < r]^ ct < 2.5, were retained. 

Leading neutron events were selected from the dijet sample by applying criteria described 
previously [20]. The main requirements are listed here. Events were required to have 
energy deposits in the FNC with energy E F ^c > 184 GeV {xl > 0.2) and timing consistent 
with the triggered event. In addition the deposits had to be close to the zero-degree point 
in order to reject protons bent into the FNC top section. Electromagnetic showers from 
photons were rejected by requiring the energy sharing among the towers to be consistent 
with a hadronic shower. Showers which started in dead material upstream of the FNC were 
rejected by requiring that the veto counter had a signal of less than one mip. Additional 
information from the FNT was used to select a subsample of events where a good position 
and thus p\ measurement was possible. The channel with the largest pulse-height in each 
of the hodoscope planes was required to be above a threshold to select neutrons which 
showered in front of the FNT plane, and transverse shower profiles were required to have 
only one peak to minimize the influence of shower fluctuations. 

After the requirements described above, the final dijet sample contained 583168 events, 
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of which a subsample of 9193 events had a neutron tag, and 4623 of these also had a well 
measured neutron position. 

The fractions of the photon and proton four-momenta entering the hard scattering, x 7 
and x p respectively, were reconstructed via 



' ~ 2E e y m 



pjetl Jetl jet2 Jet2 

Xp - 2E P ' (5j 

where r/J etl ( 2 ) and E^ n ^ are the pseudorapidity and transverse energy, respectively, of 
the highest (second highest) jet. The observable i° BS was used to separate the 
underlying photon processes since it is small (large) for resolved (direct) processes. The 
fraction of the exchanged pion four-momentum entering the hard scattering, x n in Fig. (TJ 
was reconstructed as i° BS = x° BS /(l — xl). 



4 Monte Carlo simulations 

4.1 Detector corrections 

Samples of MC events were generated to study the response of the central detector to 
jets of hadrons and the response of the forward neutron detectors. The acceptances of 
the central and forward detectors are independent and the overall acceptance factorizes 
as the product of the two; they were evaluated using two separate MC programs. 

The programs Pythia 6.221 [44] and Herwig 6.1 [45] were used to generate photopro- 
duction events for resolved and direct processes producing dijets in the central detector. 
Fragmentation into hadrons was performed using the Lund string model [46] as imple- 
mented in Jetset [47, 48] in the case of Pythia, and a cluster model [49] in the case of 
Herwig. The generated events were passed through the Geant 3.13-based [50] ZEUS 
detector- and trigger-simulation programs [29] . They were reconstructed and analyzed by 
the same program chain as the data. 

The Pythia program was used to determine the central-detector acceptance corrections. 
Samples of resolved and direct processes were generated separately. The resolved sample 
was reweighted as a function of x 7 and the direct sample as a function of W . The 
reweighting and relative contributions of the two samples were adjusted to give the best 
description of the measured x 7 and W distributions. Different reweighting and mixing 
factors were applied for the inclusive and neutron-tagged jet samples. 
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The Herwig program was used to check the systematic effects of the detector corrections. 
Direct and resolved photon processes were generated with default parameters and multiple 
interactions turned on. 

A detailed description of the efficiencies and correction factors for the leading neutron 
measurements is given elsewhere [20]. 

4.2 Model comparisons 

Previous studies have shown that MC models generating leading neutrons from the frag- 
mentation of the proton remnant do not describe the neutron xl and p\ distributions 
in DIS nor in photoproduction [20]. Models incorporating pion exchange gave the best 
description of the leading neutrons; also models with soft color interactions (SCI) [51] 
were superior to the fragmentation models. Monte Carlo programs incorporating these 
non-perturbative processes were used for comparison to the present dijet photoproduction 
data. 

The Rapgap model incorporates pion exchange to simulate leading baryon production. 
It also includes Pomeron exchange to simulate diffractive events. These processes are 
mixed with standard fragmentation according to their respective cross sections. The 
PDF parameterizations used were CTEQ5L [52] for the proton, the GRV-G LO [53] for 
the photon, the HI fit 5 [54] for the Pomeron and GRV-P LO fit [55] for the pion. The 
light-cone exponential flux factor [56] was used to model pion exchange. 

The SCI model assumes that soft color exchanges give variations in the topology of the 
confining color-string fields which then hadronize into a final state which can include a 
leading neutron. It was interfaced to the Pythia program [57]; this implementation of 
Pythia did not include multiple parton interactions. 

5 Systematic uncertainties 

Systematic uncertainties associated with the CTD and the CAL influence the jet mea- 
surement; those associated with the FNC influence the neutron measurement. They are 
considered separately. 

For the jet measurements, the systematic effects are grouped into the following classes, 
their contributions to the uncertainties on the cross sections being given in parentheses: 

• knowledge of absolute jet energy scale to 3%: (1-6%); 

• model dependence: the acceptances were estimated using Herwig instead of Pythia 
tuned as described in the previous section (5-9%); 
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• event selection: variation of W and E^ cuts by one standard deviation of the resolu- 
tion (1-6% each for W and ). 

Together, these effects resulted in uncertainties of 7-15% on the jet cross sections. The 
overall normalization has an additional uncertainty of 2.25% due to the uncertainty in 
the luminosity measurement. 

An extensive discussion of the systematic effects related to the neutron measurement is 
given elsewhere [20] ; the effects are summarized here. The neutron acceptance is affected 
by uncertainties in the beam zero-degree point and the dead material map, and uncertain- 
ties in the p\ distributions which enter into the computation of the neutron acceptance. 
The 2% uncertainty on the FNC energy scale also affects the xl and j)\ distributions. 
Systematic uncertainties from these effects were typically 5-10% of the measured quan- 
tities, for example the exponential "p\ slopes. The systematic variations largely affect 
the neutron acceptance and result in a correlated shift of neutron yields. Corrections for 
efficiency of the cuts and backgrounds in the leading neutron sample were applied to the 
normalization of the neutron yields. The corrections accounted for veto counter over- and 
under-efficiency and neutrons from proton beam-gas interactions. The overall systematic 
uncertainty on the normalization of the neutron cross sections from these corrections was 
±2.1%. Combined with the other neutron systematics, the overall systematic uncertainty 
on the total neutron rate was ±3%. 



6 Results 

6.1 Jet cross sections and ratios 

The inclusive dijet and neutron-tagged dijet phot opro duct ion cross sections have been 
measured for jets with E^ 11 ^ > 7.5(6.5) GeV and —1.5 < rf et < 2.5, in the kinematic 
region Q 2 < 1 GeV 2 and 130 < W < 280 GeV, with the additional restriction of x L > 0.2 
and 8 n < 0.75 mrad for the neutron-tagged sample. The fraction of dijet events with a 
leading neutron, the yield tln ; in the measured kinematic region is 

r LN = ae P^ e ^ Xn = 6.63 ± 0.07 (stat.) ± 0.20 (syst.)%. (6) 

In this ratio, most of the systematic effects of the dijet selection cancel, and the uncertainty 
is dominated by the systematic effects of the neutron selection. 

The differential cross sections for neutron-tagged and untagged events as functions of 
the jet variables E^ and rf et are presented in Fig. [2] and summarized in Table [U They 
contain two entries per event, one for each jet. Also shown are the neutron yields t-ln as 
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defined in Eq. (EJ as a function of the relevant variable. The cross sections as functions 
of show a reduction of about three orders of magnitude within the measured range. 
The neutron yield is approximately constant as a function of . The cross sections as 
functions of rf et rise over the range —1.5 < rf et < 0.5; for higher values of rf et they flatten. 
The neutron yield decreases with rf et . 

Figure [2] also shows the predictions of the Rapgap and SCI programs implemented as 
described in Section 14.21 Both are close in magnitude to the inclusive data. They both 
describe the steep drop with E^ and the shape of the rf et distributions. For neutron- 
tagged events Rapgap slightly overestimates and the SCI model clearly underestimates 
the cross section. They underestimate the decrease of the neutron yield with rf et . 

The differential cross sections as functions of the event variables x° BS , W and x° BS are 
presented in Fig. [3J and summarized in Table [2J The cross sections as functions of z° BS 
show two peaks at z° BS « 0.2 and x° BS ~ 0.8 which can be attributed to the resolved- and 
direct-photon contributions, respectively. The neutron-tagged sample has a significantly 
smaller resolved contribution at low x° BS . This is seen clearly in the yield which rises by 
a factor of two from low to high z° BS . The cross sections are roughly flat as a function of 
W; the yield exhibits a mild decrease with increasing W. The measured range of x° BS is 
0.04 to 0.25 and the cross section peaks close to i° BS = 0.05. The neutron yield decreases 
by a factor of two across the range measured. 

Also shown in Fig. [3J are the predictions of the Rapgap and SCI models. Rapgap does 
not have a two-peaked structure as a function of x° BS , whereas the SCI model predicts 
the drop in cross section at central values of i° BS exhibited by the data. For the neutron- 
tagged sample, Rapgap overestimates the cross section in the resolved regime while SCI 
underestimates the cross section in the direct regime. Both models predict the relatively 
weak dependence of the cross section on W and describe reasonably well the shape of 
the x° BS distribution. Neither model can reproduce the dependence of the neutron yield 
on x° BS and W. The Rapgap model predicts a small decrease of the neutron yield with 
Xp BS . However, the decrease is more pronounced in the data. The SCI model does not 
reproduce this feature at all. 

The dependence of the neutron yield on rf ct , x° BS and i° BS as seen in Figs. [2] and [3J 
indicates a violation of vertex factorization. This might be explained by the Regge factor- 
ization as discussed in Section [TJ The factorization violations seen in different variables 
are connected. A strong anticorrelation between the direct contribution (x° BS > 0.75) and 
rf ct and x° BS is apparent in the data in Fig. HI Events with low values of these variables 
contain up to 80% direct component, events with high values contain up to 90% resolved 
component. The observed drop of neutron yields at high rf et and x° BS can thus be ac- 
counted for by a lower neutron yield in the resolved photon contribution. The smaller 
dependence of the neutron yield on E 1 ^ 1 and W is consistent with this mechanism. 
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The HI collaboration has also reported similar measurements [18]. They were made in a 
similar region of , 77 jet and W as the present analysis, but restricted to xl > 0.61. The 
same pattern of vertex factorization violation was observed there. Also, after accounting 
for the different xl ranges, the cross sections are consistent. 



6.2 Neutron xl distribution and pion structure 

Figure [5] shows the normalized differential cross-section (l/a ep ^ e jjx)da ep ^> e jj X n/dxL for 
neutrons with 9 n < 0.75 mrad, which corresponds to p\ < 0.476 x\ GeV 2 . The distribution 
rises from the lowest Xl values due to the increase in p\ phase space. It reaches a 
maximum for xl ~ 0.6, and falls to zero at the endpoint xj_, — 1. Also shown are the 
predictions of the MC models. The Rapgap program gives a fair description of both 
the shape and normalization of the data, although its prediction is significantly above 
the data for xl < 0.7. The SCI model does not describe the data, predicting too few 
events with neutrons and with a spectrum peaked at too low xl. Also shown in Fig. [5] is 
the pion-exchange contribution to the Rapgap prediction for the xl distribution. This 
contribution is essential for the Rapgap prediction to describe the measured distribution. 
It dominates for xl > 0.6. Thus, in this region the dijet photoproduction data are sensitive 
to the pion structure. 

Figure [6] shows the neutron cross section as a function of log 10 (x° BS ) for xl > 0.6; the 
values are listed in Table [3j The range in i° BS is from 0.01 to 0.6; the distribution peaks 
near i° BS 0.13. Also shown in Fig. [6] are the predictions of Rapgap and SCI. The 
former provides a good description of the data while the latter underestimates the cross 
section by about a factor of three. It should be noted that Rapgap, using the pion PDF 
parameterization GRV-P LO [55] based on fixed-target data with x n > 0.1, is able to 
describe the cross section down to x^ ~ 0.01. 



6.3 p\ distributions 

The p\ distributions of the leading neutrons in different x^ bins are shown in Fig. [7| and 
summarized in Table HI They are presented as normalized doubly differential distributions, 
(1/ 'cr ep ^ e jjx)d 2 (j ep _ >e jjXn/ ' dx^dp\. The bins in p\ are at least as large as the resolution, 
which is dominated by the pr spread of the proton beam. The varying p\ ranges of the 
data are due to the aperture limitation. The line on each plot is a fit to the functional form 
da ep -> e jjXn/ ' dp\ oc exp(— hp\\ Each distribution is compatible with a single exponential 
within the statistical uncertainties. Thus, with the parameterization 

1 d 2 a ep ^ ejjXn = , x -b[x L )p 2 T 



> ejj x dx L dp\ 



a(x L ) e-^VPr, ( 7 ) 
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the neutron {xliVt) distribution is characterized by the slopes b(xi) and intercepts 
Q>{xl) = (^/o'ep~->ejjx)d 2 a e p^ e jjXn/dxLdp'^\ p 2 =Q . The results of exponential fits in bins 
of xl for the intercepts and the slopes are shown in Fig. [8] and summarized in Table [5j 
The systematic uncertainties were evaluated by making the variations discussed in Sec- 
tion [5] and repeating the fits. The intercepts fall rapidly from the lowest xl, drop mildly 
in the region xl = 0.5 — 0.8, and fall to zero at the endpoint xl — 1. In the lowest xl 
bin, the slope is consistent with zero and is not plotted; above xl = 0.5 the slope rises 
roughly linearly to a value of b ~ 13 GeV -2 at xl = 0.93. 

6.4 Comparisons of different processes 

6.4.1 Comparison to neutron production in DIS 

Figure [9] shows the normalized xl distribution of leading neutrons in dijet photopro- 
duction and in inclusive DIS with Q 2 > 2 GeV 2 [20]. The yield of neutrons from dijet 
photoproduction agrees with that in DIS at low xl < 0.4, but is lower at higher xl- For 
Xl > 0.8 the yield in dijet photoproduction is more than a factor of two lower than in 
inclusive DIS. 

Figure [9] also shows the predictions of Rapgap for dijet photoproduction and DIS. The 
predicted shapes are in fair agreement with the measurements. However, the predicted 
neutron yield is ~ 10% too high for dijet photoproduction and pa 30% too high for DIS. 
The shapes of the distributions for the two processes are compared using the ratio 

^ (1/ ^ep—^ejjXjdo^ep—^ejjXn/ dx^ .g. 

(l/cr ep ^ eX )da ep ^ eX n/dx L (Q 2 > 2 GeV 2 ) 

The result is shown in Fig. [TUJ After normalizing each prediction to its respective data 
set, Rapgap provides a fair description of the drop of the neutron yield with xl in dijet 
photoproduction relative to that in DIS. 

Figure [TT1 shows the exponential p\ slopes b(xi) for dijet photoproduction and inclusive 
DIS. They are similar in magnitude and both rise with xl- Although the slopes rise some- 
what faster with xl in the dijet photoproduction data, there is no statistically significant 
difference between the two sets except for xl > 0.9. 

6.4.2 Comparison of dijet direct and resolved photon contribu- 
tions 

The neutron xl distributions in the dijet photoproduction data, enriched in direct (x° BS > 
0.75) and resolved (x° BS < 0.75) processes, are shown in Fig. [T2J normalized to their corre- 
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sponding samples without a neutron requirement. In the resolved contribution, relatively 
fewer neutrons are observed. Figure [12] also shows the predictions of Rapgap for the 
xl distributions of the direct and resolved contributions. Figure [13] presents the ratio 
between the resolved and direct contributions to the cross section, 

_ (l/cr ep _> ejjX )da ep ^ ejj xn/dx L (x° BS < 0.75) 
PR/D (l/a ep -> ejj x)da ep -> ejj xn/dx L (x° BS > 0.75) ' 

as a function of xl for data and the Rapgap prediction. The magnitude and shape are 
not described by Rapgap. 



6.5 Role of kinematic constraints 

The xl distributions for dijet photoproduction and for DIS are depicted in Fig. O It 
is interesting to investigate whether the difference between the two distributions is a 
characteristic of the p — > n transition or if it is a kinematic effect, due to different forward 
energy flows. To investigate such kinematic constraints, X BP , the fraction of the proton 
beam energy going into the forward beampipe region, 77 > 5, was considered: 

X BP =1- (10) 

Here E p = 920 GeV is the proton beam energy and E + Pz is the longitudinal energy- 
momentum, E + Pz = ^2iEi(l + cos#i), with the sum running over all CAL and FPC 
cells with energy E{ and polar angle 6%. The energy of the leading neutron in an event is 
restricted to xl < X BP . 

The X BP distributions for dijet photoproduction and DIS, both without a leading neu- 
tron requirement, are shown in Fig. [141 The dijet photoproduction data are peaked at 
significantly lower Xbp and have a much larger tail at very low X BP than the DIS data. 
Figure [TBI shows the neutron xl distribution^ of the dijet photoproduction and DIS data 
in bins of X BP , normalized by the number of events without a neutron requirement in 
the X BP bin. They reflect the constraint xl < X BP . For any given value of X BP , the two 
samples have nearly identical xl distributions, both in shape and normalization. This 
indicates that a given value of longitudinal energy-momentum measured in the central 
detector is associated with the same neutron yield and spectrum, regardless of whether 
the process is dijet photoproduction or DIS. 



2 These xl distributions are not corrected for acceptance. The acceptance correction at a given xl 
depends only on the exponential p^, slope b{xi). As shown in Fig. [TT] the slopes for dijet photopro- 
duction and DIS have very similar values. Differences in the acceptance correction are small and may 
be ignored for the comparisons made here. 
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The effect of kinematic constraints from energy distributions in the central detector can 
also be investigated in the xl distributions of direct and resolved photoproduction as 
shown in Fig. [121 Figure [16] shows the Xbp distributions for the contributions from direct 
and resolved photons without a neutron tag being required. The resolved contribution 
peaks at lower Xbp and has a much larger tail at very low Xbp than the direct contribu- 
tion. Figure d7| shows the neutron yield as a function of xl (not corrected for acceptance) 
in different bins of Xbp for the two contributions. As in the comparison to DIS, they 
verify the constraint xl < Xbp, and for any given value of Xbp, the two samples have 
nearly identical Xl distributions, both in shape and normalization. Thus the neutron Xl 
spectra in dijet photoproduction as well as in DIS seem to depend only on the energy 
available in the proton-remnant region. 



7 Discussion of rescattering 

The good statistical accuracy of the data allows an investigation into effects of rescatter- 
ing. The comparison of photoproduction to DIS offers one way to investigate rescattering 
effects, which are predicted to result in a lower neutron yield in photoproduction. Fig- 
ure [18] shows the neutron yield as a function of xl for dijet photoproduction, inclusive 
DIS with Q 2 > 2GeV 2 , and inclusive photoproduction ep — > eXn [20]. The inclusive 
photoproduction sample was obtained by tagging the scattered positron, with a resulting 
range of Q 2 < 0.02 GeV 2 . The neutron yield for the positron-tagged inclusive sample 
agrees with the yield observed for inclusive DIS at high values of xl. At low xl, the neu- 
tron yield in inclusive photoproduction is smaller than in inclusive DIS. This was shown 
to be consistent with models of rescattering [22-25]. The neutron yield is also smaller in 
dijet photoproduction, but the xl dependence of the suppression is reversed. The neutron 
yields are similar at low values of xl, whereas the neutron yield in dijet photoproduction 
is lower at high xl values. This was shown in Fig. [TU1 

The behavior of the neutron yield for dijet photoproduction is inconsistent with the rescat- 
tering models that described the yield for the positron-tagged photoproduction sample. 
Information concerning rescattering might be difficult to obtain from a direct comparison 
of dijet photoproduction and inclusive DIS data because of the different hadronic final 
states. A Regge factorization model without rescattering effects (Rapgap) can reproduce 
reasonably the differences in neutron yields between dijet photoproduction and inclusive 
DIS. A qualitative explanation can be deduced from Eq. (0Q): The cross section is pro- 
portional to <T en _ e x(s'), and this cross section rises steeply with s' = (1 — xl)s for dijet 
production [43], whereas the cross section for the inclusive DIS reaction depends only 
weakly on s' [58,59]. Therefore one expects a drop of the ratio of the dijet to DIS neutron 
yields as s' oc (1 — xl) goes to 0. This is seen in Fig. [TUJ 
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Another way to look for such effects is to compare the direct and resolved contributions 
to dijet photoproduction. For the direct photon contribution, the photon is assumed to 
be pointlike; for the resolved photon interactions, the photon is assumed to have size and 
structure. This structure may be expected to increase the probability of rescattering. 
Indeed, the lower neutron yield in the resolved contribution to the cross section, as shown 
in the x° BS distribution in Fig. [3j seems to indicate such a loss mechanism. However, this 
seems in contradiction with the xl dependence of the effect, as shown in Figs. [T2land[T3l 
These figures show that the neutron yield in the resolved contribution decreases relative 
to the yield in the direct contribution for increasing values of Xl- This contradicts the 
predictions from the rescattering models which described the behavior of the inclusive 
photoproduction sample [20], where the effect goes in the opposite direction. Again, a 
comparison is complicated by the different hadronic final states in the direct and resolved 
contributions. 

In summary, no clear conclusion on the presence of rescattering effects in dijet photopro- 
duction can be drawn from the data alone. Only a comparison to a specific model could 
clarify this issue. 



8 Summary 

Differential cross sections for neutron-tagged and untagged dijet photoproduction, e + + 
p — > e + + jet + jet + X(+ri), have been measured. The measurements required jets 
with E^ 1 > 7.5 GeV, E% t2 > 6.5 GeV and -1.5 < rf et < 2.5, in the kinematic region 
Q 2 < 1 GeV 2 and 130 < W < 280 GeV, with the additional restriction of x L > 0.2 and 
9 n < 0.75 mrad on the neutron-tagged sample. The cross sections were measured as func- 
tions of E%\ r/J et , W, x° BS and i° BS . 

The ratios of the neutron-tagged to untagged differential cross sections show a reduction 
of the neutron yield at low x° BS , large 77 jet , and large i° BS . These regions are dominated 
by resolved photon events. 

The normalized leading- neutron xl distribution was measured. It is in reasonable agree- 
ment with the Rapgap MC model including pion exchange, which is essential to obtain 
a reasonable description of the data. In addition, the leading-neutron cross section as a 
function of x° BS , the fraction of the exchanged pion four- momentum entering the hard 
scattering, was measured in the restricted kinematic range Xl > 0.6, where pion exchange 
is the dominant production process, and good agreement with the model was found. 

The leading-neutron cross sections as a function of p\ in different regions of xl were 
measured in dijet photoproduction. The p\ distributions are well described by exponen- 
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tials, and the two-dimensional (xl,Pt) distribution is fully characterized by the slopes 
and intercepts from exponential fits in each xl bin. 

The relation between the neutron yield and the fraction of the proton beam energy going 
into the forward beam pipe region, Xbp, was studied. The relative neutron rate as a 
function of xl seems to depend only on X BP . This effect accounts for the observed 
differences between the xl distributions of the photoproduction dijet and the DIS data 
samples, and between those of the direct and resolved dijet samples. 

No clear conclusion on the presence of rescattering effects can be drawn. While the 
reduction of the neutron yield in the region enriched in resolved photons is suggestive of 
the presence of a rescattering effect, the fact that this yield reduction is mainly at large 
xl seems to contradict the basic expectations of rescattering models. 
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n 09^,^ ± 06+ 0015 + - 0029 
U.UZOO =C U.UUUO_ 0020 _ .0033 


o.572 ± o.ooTiffi ±8:811 
0.368 ±o.oo5l8:^18:8?l 
o.i62 ± 0.00318:881 18:8?^ 
0.0737 ±o.ooi8i8:88?? 18:8m 
0.0383 ± 0.001318:88^ 18S 

01 S7 ± 00 06+ 00003 + - 0018 

U.U1CW =C U.UUUD_ oo3 -0.0017 

0046 ± 0003 +a0002 +a0007 

U.UUOU ZTZ U.UUUO_ 0002 „ .0005 

001 5 ± 00 09+ 00001 + ' 0002 

U.UUIO It U.UUUZ_o.0002 -0.0002 


6.8o±o.o9i8i818:^ 
6.95 ±o.iol8;l8 ±8:81 

a no i n i O+0.19 +0.09 
~ u - iz -0.16 -0.09 

a qq + n 17+0.20 +0.12 

U.JJ IU.ll _ .20 -0.17 

7 38 + 2 c i +a24 +a12 

' - 00 31 u - ZJ -0.24 -0.14 

7 49 ± 29+ - 17 +a06 

' ^ U - ZJ -0.17 -0.11 

6 69 ± 48+ - 29 +a10 

U.UJ ZTZ U.^±0_o.29 -0.05 

e 7R _|_ n 7Q+0.20 +0.05 
O./D =C U./ O_ .20 -0.05 



^jet 


da/drf ct (nb) 




da LN 


/drji et (nb) 


r LN (%) 


-1 


33 


0.80 ±o.oii8;^ 18:1? 





078 ± 





nnfi+o-oi4 +0.010 
uuu -0.013 -0.011 


9.69 ± 





70+I.66 +0.18 
' °-0.67 -0.21 


-1 


00 


3.o7±o.o3i8: 27 *18:f 7 





259 ± 





ni 9+O.OI6 +0.012 
Uiz -0.011 -0.012 


8.42 ± 





a-i +0.42 +0.48 
^-■--0.42 -0.98 


-0 


67 


6.io±o.o4i8: 4 H8:™ 





478 ± 





rv-i 7+0.031 +0.070 
Ui ' -0.028 -0.024 


7.84 ± 





9O+0.41 +0.52 
zo -0.13 -0.52 


-0 


33 


8.84 ±0.0518;™ ±8;S 





744 ± 





nlQ +0.068 +0.044 
uiJ -0.067 -0.057 


8.41 ± 





99+0.19 +0.32 
zz -0.19 -0.82 





00 


11.00 ±o.o5i8;118;S 





847 ± 





nlQ +0.047 +0.057 
uiy -0.044 -0.042 


7.70 ± 





1 O+0.15 +0.45 
io -0.15 -0.31 





33 


12.43 ±o.o5i81? ±8-2 





926 ± 





n 9 n +0.036 +0.071 
uzu -0.036 -0.054 


7.45 ± 





1 r+0.17 +0.49 
iu -0.17 -0.14 





67 


12.98 ±0.06183 185 





897 ± 





nlQ +0.031 +0.075 
uiy -0.040 -0.062 


6.91 ± 





1 cr+0.15 +0.44 
±o -0.21 -0.10 


1 


00 


12.46 ±0.0618;! 185 





820 ± 





n 9 n +0.035 +0.046 
uzu -0.035 -0.067 


6.58 ± 





1 R +0.13 +0.19 
iu -0.18 -0.14 


1 


33 


1184 + 06 +0 ' 87 +0 - 85 

1±.0<± ZTZ U.UU_ .86 -0.80 





779 ± 





nl Q+0.053 +0.076 
uiJ -0.053 -0.070 


6.58 ± 





i fi+0.12 +0.26 
1U -0.19 -0.32 


1 


67 


12.44 ±0.06llf ±\™ 





735 ± 





m 7+O.O8O +0.085 
Ui ' -0.081 -0.066 


5.91 ± 





1/1+0.13 +0.29 
iq: -0.14 -0.11 


2 


00 


13.53 ± 0.06ll;l tlH 





800 ± 





mQ +0.108 +0.112 
uiy -0.107 -0.098 


5.91 ± 





1/1+0.07 +0.38 
lq: -0.07 -0.19 


2 


33 


12.76 ±0.061};^ 115 





707 ± 





nl q+0.067 +0.130 
Uio -0.064 -0.081 


5.54 ± 





1/1+0.04 +0.62 
± ^-0.04 -0.26 



Table 1: Differential cross-sections a (ln) for the processes e + +p — > e + ±jet±jet± 
X(+n) and the ratio tln = o"ln/o" as functions of Et andrj. For each cross section 
and ratio, the first uncertainty is statistical, the second systematic, excluding the 
CAL energy scale, and the third the systematic due to the CAL energy scale. 
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da/dx° BS (nb) 



da^/dx° BS (nb) 



r LN (%) 



0.07 
0.21 
0.36 
0.50 
0.64 
0.79 
0.93 



15.54 ± 0.11±| 
23.22 dz 0.12±| 
17.13 ±0.10lJ 
14.92 ±0.091} 
17.09 ±0.10_ s 
28.39 ±0.13+} 
21.35 ±0.11+2 



84 +1.49 
48 -2.07 
56 +1.88 
28 -2.20 
60 +1.18 
56 -1.24 

31 +1.03 

32 -1.04 

2.18 +1.17 



18 
76 
77 
16 
07 



1.21 
-1.60 
-1.84 
-0.61 
-1.00 



0.586 ±0.024 
1.174 ±0.034 
1.017 ±0.031 
1.060 ±0.032 



1.283 ± 0.037±g 
2.317 ±0.052+° 
1.949 ± 0.046+° 



103 +0.122 
088 -0.111 

137 +0.162 
110 -0.148 
088 +0.102 
086 -0.068 
091 +0.052 
098 -0.079 

164 +0.075 

165 -0.084 

151 +0.148 

138 -0.174 
158 +0.069 
132 -0.085 



3.77 
5.06 
5.94 
7.10 
7.51 
8.16 
9.13 



+ 16+ 012 
±0.15+^ 
±0.181^ 
± 0.22+°;* 
± 0.22lg;l 

+0.04 
-0.04 



±0.19 



— (— n 22+ - 34 



+0.51 
-0.29 
+0.41 
-0.17 
+0.39 
-0.19 
+0.30 
-0.54 

+0.56 
-0.32 

+0.36 
-0.45 
+0.24 
-0.11 



W (GeV) 


da/dW (n 


142 


0.109 ±0.001 


167 


0.137 ±0.001 


192 


0.143 ±0.001 


217 


0.140 ±0.001 


242 


0.132 ±0.001 


267 


0.127 ±0.001 



da LN /dW (nb/GeV) 



r-LN 7o 



+0.008 
-0.008 
+0.011 
-0.011 
+0.009 
-0.009 
+0.013 
-0.013 
+0.007 
-0.007 
+0.007 
-0.007 



+0.009 
-0.005 
+0.012 
-0.010 
+0.011 
-0.011 
+0.010 
-0.009 
+0.012 
-0.013 
+0.014 
-0.016 



0.0090 ± 
0.0097 ± 
0.0095 ± 
0.0089 ± 
0.0088 ± 
0.0078 ± 



0.0003 
0.0002 



-0.0006 
-0.0007 
-0.0008 
-0.0008 



02+ - 0006 
02+ - 0008 



0.0002 
0.0002 



-0.0005 
-0.0005 
-0.0004 
-0.0004 



-0.0012 
-0.0006 
-0.0013 
-0.0009 

-0.0011 
-0.0009 
-0.0009 
-0.0007 
-0.0007 
-0.0007 
-0.0007 
-0.0007 



8.30 ±0 
7.10 ±0 
6.62 ±0 
6.35 ±0 
6.66 ±0 
6.16 ±0 



.241°; 



18 
17 



17+° 
1 ' -o 



17 
17 



07 +0.43 

07 -0.25 

08 +0.28 
07 -0.12 

01 +0.37 
01 -0.17 
14 +0.36 
14 -0.29 
21 +0.21 
21 -0.06 
12 +0.31 
12 -0.32 



logiofa. 



OBS^ 
V > 



da/dlog l0 (x° BS ) (nb) 



da LN /dlog 10 (a;° BS ) (nb) 



r LN (%) 



-2.3 
-2.1 
-1.9 
-1.7 
-1.5 
-1.3 
-1.1 
-0.9 
-0.7 



I. 47 ± 
5.33 ± 

II. 29 ± 
16.59 ± 
21.43 ± 
24.51 ± 
14.43 ± 
3.16 ± 
0.38 ± 



n n9 +0.15 +0.62 
U.UZ_ 25 _ 39 

n n4+ -25 +1-24 
u - u ^-0.42 -0.66 

06+°- 51 + 77 

u - uu -0.57 -0.68 

n n«+°- 45 + - 99 

U.UO_ 46 _ 1Q2 

n i n +1.88 +1.75 
u - iu -1.92 -1.52 
n -I i +2.27 +2.92 
u - ii -2.23 -2.13 

n no+ 1 - 06 + 1 - 62 

U - UJ -1.02 -1.87 

n rU+°- 29 +°- 35 
u - u ^-0.29 -0.35 

n m +0.03 +0.05 
u - ul -0.03 -0.05 



0.154 ± 

0.488 ± 
0.861 ± 
1.262 ± 
1.475 ± 
1.558 ± 
0.766 ± 
0.145 ± 
0.018 ± 



0.009 



+0.014 
-0.018 



016+ ' 017 

U.U1U_ 023 

n23+ ao49 

U.UZO_ 039 

030+ - 034 

U.UOU_ 34 
u - uo ^-0.129 

n^7 +0 - 136 

u - uo ' -0.136 

n27+ a053 

u - uz ' -0.055 

n oi 1+ 016 

U.U11_q Q14 

+0.002 
-0.002 



0.004 



+0.065 
-0.023 
+0.077 
-0.041 
+0.065 
-0.056 
+0.109 
-0.069 
+0.231 
-0.161 
+0.258 
-0.171 
+0.074 
-0.090 
+0.017 
-0.015 
+0.001 
-0.001 



10.49 ± 0.62±i:S ±£S 
9.16 ± 0.3ltg %±° f 7 

7 fi o i n oi +0-17 +0.65 
1 .06 zSz U.Z1_ 17 _ 54 

7.6i ± o.i9i8:2 ±8-f2 
e.ss ± o.i6i8;| ±8:1? 



-0.19 +0.33 

-0.20 -0.21 

-0.06 +0.16 

-0.06 -0.09 

4 6 1 _l Q 37+0.21 +0.14 
^.ui =n u.o i _ 09 _ 011 



6.35 ±0.15; 
5.31 ±0.19" 



4.87 ± 1.09 



+0.35 +0.14 
-0.35 -0.77 



Table 2: Differential cross-sections o"(ln) for the processes e 



V 



■jet 



jet±X(±n) and the ratio r^ 
as in TableUi 



cln/c as functions of x° BS , W andr\. Details are 
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log 10 (x° BS ) 


da LN /dlog 10 (x° BS )(nb) 


-2 


1 


n nm n + n ooo4 +0 0005 + - 0008 

U.UU1U 3= U.UUU^±_ 0008 _ .0008 


-1 


9 


0.0226 ± 0.00321°;™! 


-1 


7 


0.111 ±0.008±g : gS?±S^g 


-1 


5 


0.248 ±0.012^°™? 


-1 


3 


432 + ni 7+0.015 +0.037 

U.^OZ 3Z U.U1 / _ .009 -0.029 


-1 


1 


^61 + 020+ 054 +0 - 070 

U.OOl 3= U.UZU_ 054 _ 33 


-0 


9 


6^)3 ~\~ 22+ 036 +0 - 082 

U.000 3= U.UZZ_ a036 _ .088 


-0 


7 


n^n + n n21 +0.060 +0.130 

U.OOU 3= U.UZ1_ 055 _ .051 


-0 


5 


0.258±0.014±g : gg±g5S? 


-0 


3 


0.0664 ±0.0068^1^ 


-0 


1 


o.oo7o ± o.oo2i±g58Sg ±g : ggg 



Table 3: Differential cross-section dai,]<j/dlog l0 (x° BS ) for the processes e + +p 
e + + jet + jet + X + n for xl > 0.6. Details are as in TableUi 
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xl range 




p 2,(GeV 2 ) 


CTnorm. (GeV~ 2 ) 


0.20-0.50 


0.38 


7.74 10" 4 
2.52 10 -3 
4.86 10" 3 
7.97 10" 3 
1.18 10~ 2 
1.65 10" 2 


1.797 ± 0.169 
1.659 ± 0.156 
1.699 ± 0.155 
1.511 ± 0.151 
1.492 ± 0.149 
1.585 ± 0.149 


0.50-0.58 


0.54 


4.84 10" 3 
1.58 10" 2 
3.03 10" 2 
4.97 10~ 2 
7.40 -10~ 2 
1.03 -10" 1 


1.135 ± 0.092 
1.008 ± 0.107 
0.808 ± 0.095 
0.915 ± 0.093 
0.884 ± 0.086 
0.694 ± 0.078 


0.58-0.66 


0.62 


6.50 -10" 3 
2.12 10" 2 
4.08 10~ 2 
6.68 10~ 2 
9.94 10 -2 
1.39 10" 1 


0.982 ± 0.073 
0.985 ± 0.091 
0.984 ± 0.089 
0.694 ± 0.070 
0.678 ± 0.064 
0.526 ± 0.058 


0.66-0.74 


0.70 


8.39 10" 3 
2.74 10" 2 
5.27 lO" 2 
8.64 lO" 2 
1.29 -10 -1 
1.79 -10" 1 


0.896 ± 0.061 
0.781 ± 0.071 
0.726 ± 0.067 
0.507 ± 0.051 
0.366 ± 0.041 
0.343 ± 0.041 


0.74-0.82 


0.78 


1.05 10- 2 
3.43 10~ 2 

6.60 -10" 2 
1.08 lO" 1 

1.61 10 _1 
2.24 -10 -1 


0.840 ± 0.053 
0.664 ± 0.058 
0.462 ± 0.047 
0.330 ± 0.036 
0.223 ± 0.028 
0.162 ± 0.024 


0.82-0.90 


0.86 


1.28 -10~ 2 
4.20 10 -2 
8.08 10~ 2 
1.33 lO" 1 
1.97 -10- 1 
2.75 10" 1 


0.364 ± 0.032 
0.289 ± 0.035 
0.194 ± 0.027 
0.145 ± 0.021 
0.044 ± 0.011 
0.048 ± 0.011 


0.90-1.00 


0.93 


1.52 KT 2 
5.03 10" 2 
9.68 -10" 2 
1.59 lO" 1 
2.36 lO" 1 
3.29 10" 1 


0.049 ± 0.009 
0.033 ± 0.009 
0.022 ± 0.007 
0.006 ± 0.003 
0.002 ± 0.002 
0.006 ± 0.003 



Table 4: The normalized doubly differential distributions o~ n 

(1/ 'a ep ^ e jjx)d 2 a ep ^ e jjXn/ ' dx^dy^. Only statistical uncertainties are shown. 
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x l ran 


ge 


(xl) 


0.20- 





50 


0.38 


0.50- 





58 


0.54 


0.58- 





66 


0.62 


0.66- 





74 


0.70 


0.74- 





82 


0.78 


0.82- 





90 


0.86 


0.90- 


1 


00 


0.93 



a ( GeV 



1.726 ± 





115 


1.084 ± 





072 


1.058 ± 





061 


0.940 ± 





054 


0.878 ± 





051 


0.420 ± 





033 


0.061 ± 





Oil 



+0.206 
-0.199 
+0.081 
-0.104 

+0.060 
-0.035 
+0.027 
-0.038 
+0.025 
-0.026 
+0.089 
-0.072 
+0.024 
-0.026 



b ( GeV -2 ) 
8.63 ± 7A5tHt 
4.00 ±1.231^ 
4.89 ± 0.83±g;|| 

+0.23 
-0.65 



6.46 ±0.72 



8.38 ± 0.67+°;^ 
9.61 ± 0.83lg;^ 
12.89 ± 2.11±|^ 



Table 5: The intercepts a and slopes b of the exponential parameterization of the 
differential cross section defined in Section ROl Statistical uncertainties are listed 
first, followed by systematic uncertainties, not including an overall normalization 
uncertainty of 2.1% on the intercepts. The systematic uncertainties are strongly 
correlated between all points. 
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Figure 1: Schematic of resolved photoproduction of dijets associated with a 
leading neutron, mediated by meson exchange. The fraction of the energy of the ex- 
changed meson (photon) participating in the partonic hard scattering that produces 
the dijet system is denoted by x n (x 7 ); the corresponding hard cross section is a. 
In direct photoproduction, the exchanged photon participates in the hard scattering 
as a point-like particle, there is no photon remnant, and x 7 = 1 . 
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Figure 2: Differential neutron-tagged and untagged dijet photoproduction cross 
sections as functions of E ] f l and rj. The ratios between cross sections, the neutron 
yields, are also given. The inner error bars, where visible, show the statistical 
uncertainty; the outer error bars, where visible, show the statistical and jet-related 
systematic uncertainties other than CAL energy scale summed in quadrature; the 
shaded bands show the contribution to the latter from the CAL energy scale. There 
is an overall systematic uncertainty on the normalization of the neutron cross- 
sections and the ratios of ±3% which is not shown. An overall uncertainty on the 
normalization of the cross sections of 2.25% due to the luminosity measurement is 
also not shown. The histograms show the predictions of the Monte Carlo models 
Rapgap (solid histogram) and Pythia with SCI (dashed histogram) as described 
in the text. 
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Figure 3: Differential neutron-tagged and untagged dijet photoproduction cross 
sections as functions of x° BS , W and log 10 (x° BS ) . The ratios between between cross 
sections, the neutron yields, are also given. Details are as in Fig. [H 
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Figure 4: Direct photon contributions (x° BS > 0.75) as functions of the other 
jet and event variables. Statistical uncertainties are smaller than the plotted solid 
points. 
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Figure 5: The normalized differential distribution (l/o~ ep ^ te jjx)do~ ep ^ te jjXn/dxL 
in dijet events. The error bars show the statistical uncertainty; neutron-related 
systematic uncertainties are shown separately as a shaded band. An overall sys- 
tematic uncertainty on the normalization of the neutron cross-sections of ±2.1% 
is not shown. The solid histogram shows the prediction of the full Rapgap model; 
the dotted histogram is the contribution from pion exchange. The dashed histogram 
is the prediction of Pythia with SCI. 
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Figure 6: Differential cross section for xl > 0.6 as a function o/log 10 (x° BS ) ; i/ie 
fraction of the exchanged pion's momentum participating in the production of the 
dijet system for the neutron-tagged sample. Details are as in Fig. The xl cut 
restricts the sample to the region where pion exchange is the dominating process. 
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Figure 7: The p\ distributions in bins of xl- The statistical uncertainties are 
shown by vertical error bars; in some cases they are smaller than the plotted symbol. 
The systematic uncertainties are not shown. The line on each plot is the result of 
a fit to the form da ep ^ e jjx n / ' dp\ oc exp(— bp%). 



29 



■r 2 
> 

O 1.6 



CO 



1.2 
0.8 
0.4 




ZEUS 



* ♦ 

ZEUS 40 pb" 1 
ep -> ejjXn 
Q 2 < 1 GeV 2 
p 2 . < 0.476 x* GeV 2 



> 
a) 

a 12 



8 



0.4 0.6 



0.8 



Systematic 
uncertainty 



♦ 



0.4 0.6 0.8 1 

x, 



Figure 8: Intercepts a and exponential slopes b versus xl from fits of the p\ 
distributions to the form (l/c ep ^ e jjx)d 2 a ep ^ e jjXn/dxLdp'^ = aexp(-bp^) over the 
range p\ < 0.476x| GeV 2 . The error bars show the statistical uncertainties; the 
shaded bands show the neutron-related systematic uncertainties. The band for the 
intercepts does not include the overall normalization uncertainty of ±2.1%. 
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Figure 9: Normalized differential distributions (1/ 'c ep ^ e x)da ep ^ e xn/ 'dx^. The 
solid points are for dijet photoproduction and the open points for DIS [20]. Both 
distributions are normalized by their respective untagged cross sections. Statistical 
uncertainties are shown as vertical bars; in the DIS case they are smaller than 
the plotted symbols. The systematic uncertainties, shown as the shaded band, are 
similar for both data sets. The histograms are the predictions o/RAPGAP. 
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Figure 10: Ratio of the leading-neutron xl distributions dijet photoproduction 
to DIS. The data are the solid points, with statistical uncertainties shown by ver- 
tical bars and the systematic uncertainty by the shaded band. The histogram is 
the prediction of Rapgap after its normalization was adjusted to both data sets 
separately. 
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Figure 11: Exponential slopes b versus xl from fits of the p\ distributions to the 
form da ep ^ e xn/ dpr ^ ex P( — fy°r) over the kinematic range p\ < 0. 476:c! GeV 2 . 
The solid points are for dijet photoproduction, the open points for DIS. Statistical 
uncertainties are shown as vertical bars, where visible. The systematic uncertain- 
ties, shown as the shaded band, are similar for both data sets. 
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Figure 12: The normalized differential distributions {l/o~ ep ^ e jjx)da ep ^ e jjXn/dxL 
for the direct- enhanced (x® > 0.75, solid points) and resolved- enhanced (x® < 
0.75, open points) dijet photoproduction samples. Statistical uncertainties are 
shown as vertical bars. The systematic uncertainties, shown as the shaded band, 
are similar for both data sets. The histograms are the predictions of Rapgap for 
the respective components. 



34 



ZEUS 



□ 1.2 

Q. 



0.8 



0.6 



0.4 



0.2 



I t i 



t t 



ZEUS 40 pb" 1 

ep -> ejjXn 

Q 2 < 1 GeV 2 

p 2 . < 0.476 x* GeV 2 

RAPGAP 



i 



I 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x, 



Figure 13: Ratio of leading-neutron spectra for resolved (x® BS < 0.75) and direct 
(x® BS > 0.75) contributions to dijet photoproduction. Only statistical uncertainties 
are displayed. The data is compared to the prediction o/Rapgap. 
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Figure 14: Comparison between Xbp distributions for the dijet photoproduc- 
tion (solid points) and inclusive DIS (open points) [20] samples. In both cases 
no neutron tag was required. Statistical uncertainties are smaller than the plotted 
points. 



36 



ZEUS 



_i 0.04 



x 

T3 



c 

^ 0.02 



Q. 



\ o 

Q) 



0.55<X BP <0.6 



1 



O.6<X BP <O.05 
9 



: • • a o — 



a • 



0.65<X BP <0.7 



O ! 



! » O tt 



T 



0.06 


0.7<X BP <0.75 


i i 


i i i i 

0.75<X BP <0.8 
* 


i 


i i i i i i i 

0.8<X BP <0.85 * j 
o : 

* 


0.04 


<5 o 









» : 

• 


0.02 










• 

:* * : 


n 


— i i 1 i i ■ i O i O- 


■ — i i 


• <^ 





0.075 
0.05 
0.025 




1 1 1 1 1 

0.85<X BP <0.9 


_ 


i i i i i i i _ 

0.9<X BP <0.95 ♦ j 


l l l l l l 

0.95<X BP <1 


a ♦ 






• dijet Q 2 < 1 GeV 2 


• <* ° 




Q A 

-a 


O DIS Q 2 > 2 GeV 2 


o 




« * * ' 


— ill 


9- 


— i i i i i i ■ — 





0.5 



0.5 



0.5 



Figure 15: Neutron yield as a function of xl for different bins of Xbp for 
the dijet photoproduction (solid points) and inclusive DIS (open points) samples. 
The data are not corrected for detector acceptance. Statistical uncertainties are 
shown as vertical bars, where visible. The vertical dashed lines show the constraint 
x L < X BP . 
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Figure 16: Comparison between Xbp distributions for the dijet photoproduction 
direct (solid points) and resolved (open points) photon contributions. In both cases 
no neutron tag was required. Statistical uncertainties are smaller than the plotted 
points. 
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Figure 17: Neutron yield as a function of Xl for different bins of Xbp for 
the dijet photoproduction direct (solid points) and resolved (open points) photon 
contributions. The data are not corrected for detector acceptance. Statistical un- 
certainties are shown as vertical bars, where visible. The vertical dashed lines show 
the constraint xl < X BP . 
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Figure 18: Neutron yields as a function of Xl for dijet photoproduction (solid 
points), inclusive DIS (open points), and inclusive photoproduction (shaded trian- 
gles) [20]. Statistical uncertainties are shown as vertical bars, where visible. The 
systematic uncertainties, shown as the shaded band, are similar for all three data 
sets. 
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